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We prove the existence of the exact density-functional theory formalism for open electronic sys- 
tems, and develop subsequently an exact time-dependent density-functional theory (TDDFT) for- 
mulation for the dynamic response. The TDDFT formulation depends in principle only on the 
electron density of the reduced system. Based on the nonequilibrium Green's function technique, 
it is expressed in the form of the equation of motion for the reduced single-electron density ma- 
trix, and this provides thus an efficient numerical approach to calculate the dynamic properties of 
open electronic systems. In the steady-state limit, the conventional first-principles nonequilibrium 
Green's function formulation for the current is recovered. 
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Density-functional theory (DFT) has been widely used 
as a research tool in condensed matter physics, chemistry, 
materials science, and nanoscience. The Hohenberg- 
Kohn theorem pj lays the foundation of DFT. The Kohn- 
Sham formalism [2( provides the practical solution to 
calculate the ground state properties of electronic sys- 
tems. Runge and Gross extended further DFT to calcu- 
late the time-dependent properties and hence the excited 
state properties of any electronic systems 0- The accu- 
racy of DFT or TDDFT is determined by the exchange- 
correlation functional. If the exact exchange-correlation 
functional were known, the Kohn-Sham formalism would 
have provided the exact ground state properties, and 
the Runge-Gross extension, TDDFT, would have yielded 
the exact properties of excited states. Despite of their 
wide range of applications, DFT and TDDFT have been 
mostly limited to closed systems. 

Fundamental progress has been made in the field of 
molecular electronics recently. DFT-based simulations 
on quantum transport through individual molecules at- 
tached to electrodes offer guidance for the design of prac- 
tical devices [_D- 13- E3| - These simulations focus on the 
steady-state currents under the bias voltages. Two types 
of approaches have been adopted. One is the Lippmann- 
Schwinger formalism by Lang and coworkers Q- The 
other is the first-principles nonequilibrium Green's func- 
tion technique j^, 0, 0, 0, [^. In both approaches 
the Kohn-Sham Fock operator is taken as the effective 
single-electron model Hamiltonian, and the transmission 
coefficients are calculated within the noninteracting elec- 
tron model. It is thus not clear whether the two ap- 
proaches are rigorous. Recently Stefanucci and Almbladh 
derived an exact expression for time-dependent current 
in the framework of TDDFT . In the steady-current 
limit, their expression leads to the conventional first- 
principles nonequilibrium Green's function formalism if 
the TDDFT exchange-correlation functional is adopted. 
However, they did not provide a feasible numerical formu- 
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lation for simulating the transient response of molecular 
electronic devices. In this communication, we present a 
rigorous first-principles formulation to calculate the dy- 
namic properties of open electronic systems. We prove 
first a theorem that the electron density distribution of 
the reduced system determines all physical properties or 
processes of the entire system. The theorem lays down 
the foundation of the first-principles method for open sys- 
tems. We present then the equation of motion (EOM) for 
nonequilibrium Green's functions (NEGF) in the frame- 
work of TDDFT. By introducing a new functional for 
the interaction between the reduced system and the en- 
vironment, we develop further a reduced-single-electron- 
density-matrix-based TDDFT formulation. Finally, we 
derive an exact expression for the current which leads 
to the existing DFT-NEGF formula in the steady-state 
limit. This shows that the conventional DFT-NEGF for- 
malism can be exact so long as the correct exchange- 
correlation functional is adopted. 

Both Hohenberg-Kohn theorem and Runge-Gross ex- 
tension apply to isolated systems. Applying Hohenberg- 
Kohn-Sham's DFT and Runge-Gross's TDDFT to open 
systems requires in principle the knowledge of the elec- 
tron density distribution of the total system which con- 
sists of the reduced system and the environment. This 
presents a major obstacle in simulating the dynamic pro- 
cesses of open systems. Our objective is to develop an 
exact DFT formulation for open systems. In fact, we 
are interested only in the physical properties and pro- 
cesses of the reduced system. The environment provides 
the boundary conditions and serves as the current source 
and energy sink. We thus concentrate on the reduced 
system. 

Any electron density distribution function p(r) of a 
real physical system is a real analytic function. We may 
treat nuclei as point charges, and this would only lead to 
non-analytic electron density at isolated points. In prac- 
tical quantum mechanical simulations, analytic functions 
such as Gaussian functions and plane wave functions are 
adopted as basis sets, which results in analytic electron 
density distribution. Therefore, we conclude that any 
electron density functions of real systems are real ana- 
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lytic on connected physical space. Based on this, we show 
below that for a real physical system the electron density 
distribution function on a sub-space determines uniquely 
its values on the entire physical space. This is nothing 
but the analytic continuation of a real analytic function. 
The proof for the univariable real analytical functions can 
be found in textbooks, for instance, reference The 
extension to the multivariable real analytical functions is 
straightforward. 

Lemma: The electron density distribution function 
p(r) is real analytic on a connected physical space U. 
W C U is a sub-space. If p(r) is known for all r 6 W, 
p(r) can be uniquely determined on entire U . 

Proof: To facilitate our discussion, the following nota- 
tions are introduced. Set Z+ = {0, 1,2,.. .}, and 7 is an 
element of (Z+) 3 , i.e., 7 = (71,72,73) S (Z+) 3 . The dis- 
placement vector r is denoted by the three-dimensional 
variable x = (x±, X2, #3) G U. Denote that 7! 

Q~<1 Q^S 
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and 
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Suppose that another density distribution function 
p'(x) is real analytic on U and equal to p(x) for all 

x G W. We have = ^0^- for all x e IF and 

7 G (Z + ) 3 . Taking a point xo at or infinitely close 
to the boundary of W , we may expand p(x) and p{x') 



around xo, 



7 g(Z+)3 7 | 



G> T p(z) 



(x - Xo) 7 



and p'(x) = E 7 e(z+)3 ^ 



(x — Xq) 7 . Assuming 



that the convergence radii for the Taylor expansions of 
p(x) and p'(x) at xo are both larger than a positive fi- 
nite real number b, we have thus p{x) = p'(x) for all x G 

7^1 / \ r 1 ^ j t d 1 p(x) d 1 p' (x) 

Db(x ) = {x : |x - x | < b\ since & H x \ = g£ . 

x xo 

Therefore, the equality /o'(x) = p(x) has been expanded 
beyond W to include D (xo). Since U is connected the 
above procedure can be repeated until p'(x) = p(x) for 
all x G U. 

We have thus proven that p can be uniquely deter- 
mined on U once it is known on W, and are ready to 
prove the following theorem. 

Theorem: Electron density function p(r) for a sub- 
system of a connected real physical system determines 
uniquely all electronic properties of the entire system. 

Proof: Assuming the physical space spanned by the 
subsystem and the real physical system are W and U, 
respectively. W is thus a sub-space of U, i.e., W C U. 
According to the above lemma, p(r) on W determines 
uniquely its values on U, i.e., p(r) of the subsystem de- 
termines p(r) of the entire system. 

Hohenberg-Kohn theorem and Runge-Gross extension 
state that the electron density distribution of a system 
determines uniquely all its electronic properties. There- 
fore, we conclude that p(r) for a subsystem determines 
all the electronic properties of the real physical system. 

The above theorem guarantees the existence of an ex- 
act DFT-type method for open systems. In principle, all 
we need to know is the electron density of the reduced 
system. The electron density distribution in the environ- 
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FIG. 1: Schematic representation of the experimental setup 
for quantum transport through a molecular device. 



ment can be obtained by the analytic continuation of the 
electron density function at or near the boundary. The 
challenge is to develop a practical first-principles method. 

Fig. ^ depicts one type of open systems, a molecular 
device. It consists of the reduced system or device region 
D and the environment, the left and right electrodes L 
and R. Taking this as an example, we develop an ex- 
act DFT formalism for the open systems. To calculate 
the properties of a molecular device, we need only the 
electron density distribution in the device region. The 
influence of the electrodes can be determined by the elec- 
tron density distribution in the device region. Within the 
TDDFT formulation, we proceed to derive the EOM for 
the lesser Green's function: 



GZ m (t,t')=i(al(t')a n (t)) 



(1) 



where a n (t) and a^it') are the Heisenberg annihilation 
and creation operators for atomic orbitals n and m in the 
reduced system D at time t and t' , respectively. Based on 
the Keldysh formalism |15| and the analytic continuation 
rules of Langreth 0] , Jauho et al. developed a NEGF 
formulation for current evaluation 1171. Based on the 
same procedure adopted in reference [Tj| , we obtain 
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+Kl(t,r)G< m (T,t'] 



(2) 



where Gf m (T,t') is the advanced Green's function @|, 
t) and r) are the self-energies on D induced 

by L and R whose expressions can be found in references 
such as or i an( A h n i (t) is the Kohn-Sham Fock ma- 
trix element. Eq. (0) is the exact TDDFT formulation 
for the open electronic systems. However, G< m (t,t') and 
Gf m (T,t') are the two-time Green's functions. It is thus 
extremely time-consuming to solve Eq. J2J numerically. 
Alternative must be sought. 

Yokojima et al. developed a dynamic mean-field the- 
ory for dissipative interacting many-electron systems |l8j . 
An EOM for the reduced single-electron density matrix 
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was derived to simulate the excitation and nonradiative 
relaxation of a molecule embedded in a thermal bath. 
This is in analogy to our case although our environment 
is actually a fermion bath instead of a boson bath. The 
major difference is that the number of electrons in the 
reduced system is conserved in reference 0] while in our 
case it is not. Note that the reduced single-electron den- 
sity matrix a is actually the lesser Green's function of 
identical time variables, 



,(*) 



^G nm (t 7 ^ ) | t'=t 



(3) 



Thus, the EOM for a can be written down readily with 
the aid of Eq. J5J , 



d 



leD 



G nm {t, ^ / 



+ ~EI/G nm (t, * ) 



a=L,R 



t'=t 

(*), (4) 



where Q a ,nm(t) on the right-hand side (RHS) is the dis- 
sipative term due to the lead a (L or R) whose expanded 
form is 



/oo 
dr 



G<(t,T)K,lrn(r,t) 



+ G r nl (t, r)E< (r, t) - E< (t, r)G? m (r, t) 



(5) 



And the current through the interfaces Sl or Sr (see 
Fig. ^| can be expressed as 

J «w = -^(E c Lw <*«(*)) 



\k£a 



E E (*. *) - *w,k a w 



feGa 2G-D 

20? j [ drtr 
+ G r D (t,T)Z<(r,t) 



-i tr 



Qa(t) 



G<(t,r)X a a (T,t) 



*(*), 
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(6) 



where (i) is the coupling matrix element between 
the atomic orbital / and the single-electron state k a in 
£ or R, G< n /t,t>) ee i(a\{t')c ka {t)) and G< k Jt,t>) ee 



i(ct (t')ai(t)), Ckjt) and c[ (t') are the annihilation and 
creation operators for k a , respectively. At first glance 
Eq. 10} is not self-closed since the Gs are to be solved. 
According to the theorem we proved earlier, all physi- 
cal quantities are explicit or implicit functionals of the 
electron density in D, po(t)- Gs and Es are thus also 
universal functionals of po{t)- Therefore, we can recast 



Eq. (0J into a formally closed form, 



1(7 



h[ PD (t)],a\ + Qa\PD(t)]- 



(7) 



a=L,R 



Neglecting the second term on the RHS of Eq. (TJ leads 
to the conventional TDDFT formulation in terms of re- 
duced single-electron density matrix 19j. The second 
term describes the dissipative processes where electrons 
enter and leave the region D. Besides the exchange- 
correlation functional, the additional universal density 
functional Q a [pD(t)} is introduced to account for the dis- 
sipative interaction between the reduced system and its 
environment. Eq. J7J is thus the TDDFT formulation in 
terms of the reduced single-electron matrix for the open 
system. In the frozen DFT approach [2(j an additional 
exchange-correlation functional term was introduced to 
account for the exchange-correlation interaction between 
the system and the environment. This additional term is 
included in h[p D (t)\ of Eq. Q). Admittedly, Q a [pD(t)} 
can be an extremely complex functional. Progressive 
approximations are needed for the practical solution of 
Eq. Q). Compared to Eq. J5J, Eq. Q) may be much more 
convenient to be solved numerically. 

To obtain the steady-state solution of Eqs. @ or fl, 
we adopt a similar strategy as that of reference |13| . 
As t,r +oo, 1*5, (i,r) = V n ,kJt)V^ ik Jr) becomes 
asymptotically time-independent, and Gs and Es rely 
simply on the difference of the two time- variables |13| . 
The expression for the steady-state current is thus as fol- 
lows, 



Jl(oo) 



T(e) 



-Jr(oo) = - j Q L ,nn(<X>) 
neD 



[fL(e)-f R (e)] T{e)de, 
2TTT] L (e)r] R (e) 

xtr\G r D (e)T R (e)G a D (e)T L (e) 



(8) 
(9) 



Here T(e) is the transmission coefficient, f a (e) is the 
Fermi distribution function, and rj a (e) — J2kea — e fc) 
is the density of states for the lead a (L or R) . Eq. © is 
exactly the Landauer formula HHH^I in the DFT-NEGF 
formalism 0, Q . The only difference is that Eq. J3J is 
derived within the TDDFT formalism in our case while 
it is evaluated within the DFT framework in the case of 
the DFT-NEGF formulation In other words, the 

DFT-NEGF formalism can be exact so long as the cor- 
rect exchange-correlation functional is used! This is not 
surprising, and is simply a consequence of that (i) DFT 
and TDDFT can yield the exact electron density and (ii) 
the current is the time derivative of the total charge. 

Just as the exchange-correlation functional, the ex- 
act functional form of Q a on density is rather difficult 
to derive. Various approximated expressions have been 
adopted for the DFT exchange-correlation functional in 
the practical implementation. Similar strategy can be 
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employed for Q a . One such scheme is the wide-band 
limit (WBL) approximation |l7|. which consists of a se- 
ries of approximations imposed on the leads: (i) their 
band-widths are assumed to be infinitely large, (ii) their 
linewidths A£(f, t) defined by Trr] a (s%)r ka (t, t) are re- 
garded as energy independent, i.e., A%(t, r) « A a (t, r) « 
A Q , and (iii) the energy shifts are taken as level indepen- 
dent, i.e., 5e%(t) w » & Q for L or _R. The physical 
essence of the transport problem is captured under these 
reasonable hypotheses [T^l- In the practical implemen- 
tation, the effects of the specific electronic structures of 
the leads can be captured by enlarging the device region 
to include enough portions of the electrodes. 

Following the WBL approximation in reference [l7j . 
we obtain that 

Q a =[P a (t)-[P a (ttf] -i{A a ,a}, (10) 

where the curly bracket on the RHS denotes the anticom- 
mutator, and by taking t = as the switch-on instant 
P a (t) can be expressed as 

P a (t) = -U^\t){ / dr- 17<+>(t) 

7T ,/n t — T 



■i[Se a t+ l x^°{t-r)] \ 

dr 1 -^—. -G^(-t) A c 



t - T 



2iA a , 



(11) 



where ^ is the chemical potential of the lead a (L or 
R) in its initial ground state, G r ^°(— r) is the retarded 
Green's function of D before the switch-on instant, and 
Us are defined as 

U {±) {t) = exp(±?: f h{r)dT± AQ< V ( 12 ) 

^ J ° a=L,R > 



Eqs. l(TU) l — ifESj l constitute the WBL formulation of the 
TDDFT-NEGF formalism. Although its explicit func- 
tional dependency is not given, Q a depends implicitly 
on p D via Eqs. (fTTTf) - jQ . 

To summarize, we have proven the existence of the ex- 
act TDDFT formalism for the open electronic systems, 
and have proposed a TDDFT-NEGF formulation to cal- 
culate the quantum transport properties of molecular 
devices. Since TDDFT results in formally exact den- 
sity distribution, the TDDFT-NEGF formulation is in 
principle an exact theory to evaluate the transient and 
steady-state currents. In particular, the TDDFT-NEGF 
expression for the steady-state current has the exact same 
form as that of the conventional DFT-NEGF formal- 
ism 0, 0, El O, 0, and this provides rigorous the- 
oretical foundation for the existing DFT-based method- 
ologies @,H,lilElE10 calculating the steady currents 
through molecular devices. 

In addition to the conventional exchange-correlation 
functional, a new density functional is introduced to ac- 
count for the dissipative interaction between the reduced 
system and the environment. In the WBL approxima- 
tion, the new functional can be expressed in a relatively 
simple form which depends implicitly on the electron den- 
sity of the reduced system. Since the basic variable in our 
formulation is the reduce single-electron density matrix, 
the linear-scaling techniques such as that of reference 0] 
can be adopted to further speed up the computation. 
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